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Chapter 1 
Introduction 


It is the intent of this report to investigate possible treatments for a two-dimensional 
(2D) edge termination based upon both physical and electromagnetic constraints. 
In the context used here, an edge termination is a leading or trailing portion of an 
aerodynamic airfoil. The shape of the airfoil is largely controlled by aerodynamic con- 
straints. However, the use of penetrable materials allows electromagnetic constraints 
to be simultaneously realized. A variety of computational methods are available to 
evaluate electromagnetic scattering from two-dimensional edge terminations. These 
techniques are continuously being improved in method, speed, and efficiency. 

A traditional numerical approach is a method of moments (MoM) technique [1]. 
This technique can be efficiently applied to perfectly conducting or homogeneous 
penetrable structures since the unknowns can be limited to the perimeter of the 
structure. However, when the structure is inhomogeneous, the entire internal region 
of the structure has to be expanded with unknowns. This expansion causes the MoM 
approach to become exhaustive in terms of the computer resources and computation 
time required for generation and solution of the resulting linear system of equations. 

An alternate numerical approach is a hybrid combination of both the finite element 
method (FEM) and the MoM [2]. This hybrid method still requires the expansion of 
the internal region, with the difference being that a sparse matrix is generated through 
the use of the FEM. The internal fields of the FEM solution are then coupled to the 
unbounded free space region through an integral equation formulation commonly 
called the boundary element method (BEM). 
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An unfortunate consequence of electromagnetic analysis for realistic structures is 
that a significant amount of input information is required. This input information 
is a discretization of the geometry under consideration, i.e., a (2D) FEM mesh. In 
addition, the design of a satisfactory edge termination is a problem that involves many 
material and structural parameters. As a computational means of searching for and 
keeping track of parameter combinations with desirable performance, an optimization 
method is used that is termed genetic algorithm (GA) [3]. This type of algorithm, 
which is based upon a self learning scheme, has many advantages over traditional 
optimization methods. Two advantages are that no derivative information is required 
and multiple solutions with different parameter combinations are possible. This GA 
technique is adapted for this investigation to provide ease of integration between the 
optimizer and an internally developed geometry mesher. 

Chapter 2 presents a description of the theory and implementation of two types of 
triangular mesh generation schemes: Delaunay triangulation and the advancing front 
method. A description of these techniques and their implementation is provided. 
Chapter 3 presents an overview of the basic terminology and concepts of GAs. A GA 
software package called GENESIS [4] is applied to two problems. The results of the 
two problems are discussed and interpreted. 

Chapter 4 contains the results of an investigation using a variety of treatments 
that can be applied to edge terminations. Among these treatments are conducting 
bulkhead orientation and shaping, variation of dielectric and/or conducting wedge an- 
gle(s), internal material loadings, placement of resistive cards, and the nulling effects 
of material rods. Chapter 5 uses selected results of Chapter 4 as a guide in search- 
ing for edge terminations that demonstrate desirable performance over the frequency 
range of interest for both transverse magnetic and transverse electric polarizations. 
The GA is also applied to the treatments of Chapter 4 to assist in improving current 
designs, as well as to explore yet unconsidered designs. All of the FEM calculations 
in Chapters 4 and 5 are either swept frequency (echo width versus frequency) or pat- 
tern (echo width versus angle) calculations. Additionally, time domain information 
pertaining to the downrange scattering characteristics of the structure is obtained 
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from swept frequency domain calculations. The document concludes with Chapter 6 
providing a summary of the edge treatments found to be most effective. 
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Chapter 2 

Automatic Two Dimensional 
Triangular Mesh Generation 

2.1 Initial Considerations 

2.1.1 Mesh Generation Topics 

Three topics encountered when considering mesh generation are: mesh generation 
in two and three dimensions; mesh generation techniques using either quadrilat- 
eral/cubic or triangular/tetrahedral element meshes; and manual and semi-automatic 
mesh generation. For the purposes of the Finite Element Method (FEM) calculations 
that are performed for this report, an automatic two dimensional triangular mesh 
generation scheme is desired. Once this type of mesh generation is chosen, two trian- 
gulation concepts are found that are the most popular: Delaunay triangulation, and 
the advancing front method. 

Delaunay triangulation is desirable since it insures that small included angles 
within each element are avoided. An efficient method of computing the Delaunay 
triangulation using Watson’s [5] algorithm is described, complete with a FORTRAN 
source code, in [6]. This triangulation method can be incorporated into a mesh gen- 
eration code, with two disadvantages: only convex polygonal regions can be properly 
triangulated and, due to this fact, concave geometries have to be broken into convex 
polygonal regions before meshing. 

The advancing front method is desirable in that it can create meshes in convex 
and concave polygonal regions, and that the method can be used to mesh geome- 
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tries containing voids. Lo [7] provides a straightforward method of implementing the 
advancing front method, as well as a general nodal generation scheme. 

2.1.2 Basic Mesh Generation Components 

A basic mesh generator creates and provides a set of points to be connected (in this 
case triangulated), as well as the nodal connectivity of the points specifying the ele- 
ments (triangles). Additionally, a mesh generator must define geometry information 
concerning mesh boundaries, imposed boundary conditions, and material geometry 
parameters. 

The sections of this chapter will address mesh generation considerations and imple- 
mentation with the last section providing illustrations of several meshed geometries. 

2.2 Desirable Mesh Generation Characteristics 

A good summary of desired characteristics of a mesh generation procedure is found 
in [8]: 

• The algorithm should be able to handle arbitrary geometries in a fully automatic 
manner and with minimum user intervention. 

• The input data should be reduced to a computerized geometric representation 
of the domain to be discretized. 

• The approach followed should provide control over the spatial variation of ele- 
ment size and shape through the domain. 

• Adaptive methods should be incorporated into the process, with the objective of 
producing the most accurate approximation of the solution for a given number 
of points. 

The first two characteristics listed here are accommodated by the implemented 
algorithm and are described in the remainder of this chapter. Spatial variation of 
elements and adaptive methods are more advanced mesh generation topics. Since 


5 



the main focus of the work presented here is not mesh generation, a constant spa- 
tial variation of elements with no adaptive methods is used. When possible, FEM 
results calculated using this mesh generation scheme are compared with previous cal- 
culations, eigenfunction solutions, and moment method solutions to guarantee the 
solution accuracy (Several reference curves can be found in [2]). Adaptive methods 
can, however, be incorporated into the algorithm in the future should the need arise. 

2.3 Mesh Generation Philosophy 

Mesh generation can be divided into two steps: nodal generation and element creation. 
However, before these steps can take place, an adequate geometry definition is needed. 
The following sections describe the information contained in the geometry file, as well 
as the algorithms used for nodal definition and element creation. 

2.3.1 Geometry Definition 

For the geometries that will be defined in Chapter 4, a complete geometry file contains: 

• The number of vertices in the geometry. 

• Each vertex’s X and Y coordinates. 

• The number of regions in the geometry. 

• The number of boundary segments for each region. 

• Each boundary segment’s type: line or arc. If the segment is an arc, the ver- 
tex number of its center and the arc’s direction of rotation must be specified 
(clockwise or counterclockwise). 

• The vertex numbers of each boundary segment’s endpoints. 

• Each boundary segment’s enforced boundary condition, which can be one of the 
following: 

1. No enforced boundary condition 
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2. Perfect Electric Conductor (PEC) - Transverse Magnetic (TM) 

3. Perfect Electric Conductor (PEC) - Transverse Electric (TE) 

4. Resistive 

If a segment is resistive, a reference number is added following the boundary 
condition specification. The reference number corresponds to an impedance 
value contained in the material file. 

2.3.2 Nodal Generation 

Nodal generation consists of three parts: 

1. Generation of nodes on region boundaries. 

2. Compression of region boundary nodes to eliminate doubly defined nodes. 

3. Generation of nodes interior to region boundaries. 

Once all nodes are adequately defined, they are triangulated. 

2.3.3 Triangulation 

Triangular mesh elements were selected over quadrilateral elements. In general, tri- 
angular meshes are used most often as they are able to mesh any polygonal region 
exactly. Once the decision is made to use triangular elements, a method of tri- 
angulation is chosen. The two most popular methods of triangulation are Delaunay 
triangulation and the advancing front method. Both methods have been implemented 
and will be presented here. The Delaunay triangulation method is implemented with 
details and source code given in [6]. Additionally, an advancing wavefront algorithm 
is implemented based upon Lo’s [7] algorithm. 

2.4 Two Dimensional Mesh Generation Implemen- 
tation 

Two complete mesh generation programs have been written in FORTRAN. The dif- 
ference between the two codes is in their method of triangulation. The logic used in 
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the remaining sections of the two codes (excluding triangulation) is equivalent, with 
only slight differences in the logic’s implementation. The equivalent sections of code 
are: 


• ASCII interface, 

• Geometry specification and 

• Nodal generation. 

Once a geometry has been meshed, the necessary information is written to a file 
in either ASCII or binary format. 

2.4.1 ASCII Interface 

The user interface to control the mesh is provided through the follow ASCII screen 
description: 

*********** input parameters ********************* 

1- SPECIFY GEOMETRY FILE NAME: strip3 

2- ENTER NODAL SPACING: 0.1000000 

3- SPECIFY OUTPUT FILE NAME (ROOT) : strip3 

4- OUTPUT FORMAT (T,F:BINARY, ASCII) T 

5- PLOT GRID: (n,s,h) h 

************************************************** 

ANY CHANGES? (INPUT ROW # OR 0=N0) : 

The required inputs are self-explanatory. The units for the nodal spacing are 
the same as the units used to describe the geometry in the input geometry file. 
Specifying the output format concatenates an extension of “.bin” or“.asc” to the end 
of the output file name. Graphical output can be obtained through one of two means. 
Specifying an “s” generates a UNIX plot file while an “h” generates a postscript file. 
An “n” produces no graphical information. 

2.4.2 Geometry Specification 

The sample geometry considered is shown in Figure 2.1. The actual geometry input 
set, which is provided below, begins with the number of vertices in the geometry on 
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the first line of the data file, followed by a corresponding list of X/Y coordinate pairs 
for each vertex. Next, the number of regions in the geometry is listed. Finally, for 
each region of the geometry, there is a corresponding set of information. The number 
of boundary segments in the region is given at the beginning of the region informa- 
tion. The remaining region information contains each line segment’s type: ’arc’ or 
’line’; each segment’s endpoints listed in terms of their corresponding vertex numbers; 
and each segment’s enforced boundary condition: ’none’ - normal, ’pectm’ or ’pecte’ 
- perfect electric conductor, TM or TE, ’res’ - resistive card. If the segment is an arc, 
the next line contains the arc’s center vertex number and direction of rotation: ’cw’ 
or ’ccw’. If the segment’s boundary condition is resistive, a complex number corre- 
sponding to the segment’s impedance (the actual impedance is stored in an external 
material file) is added on the line following the line or arc segment information. Each 
region’s boundary segments should be defined in a counterclockwise manner. 
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It should be noted that the code that utilizes Delaunay triangulation was written 
to triangulate convex polygonal regions, with region boundaries made up only of 
straight line segments. Therefore the ’arc’/’line’ distinction as well as the arc center 
and arc rotation direction information are not included in the geometry file when 
using the Delaunay triangulation code. 


5 

0 25 
0 20 
0 0 
0 -20 
0 -25 
2 
4 

'line’ 2 1 ’none' 
’arc* 1 5 'none’ 

3 ’ccw’ 

’line’ 5 4 ’none’ 
’arc’ 4 2 ’res’ 

3 ’cw’ 

(-3,0) 

4 

’line’ 1 2 ’none’ 
’arc’ 2 4 ’pectm’ 
3 ’ cw’ 

’line’ 4 5 ’none’ 
’arc’ 5 1 ’none’ 

3 ’ccw’ 


! 5 vertices 
! vertex coordinates 


2 regions 
Reg 1, 4 segments 
straight seg, no BC 
arc seg, no BC 

center point, rotation direction 
straight seg, no BC 
arc seg, resistive 
center point, rotation direction 
material file index number 
Reg 2, 4 segments 
straight seg, no BC 
arc seg, PEC TM 

center point, rotation direction 
straight seg, no BC 
arc seg, no BC 

center point, rotation direction 


2.4.3 Nodal Generation 

Nodal Creation on Domain Boundaries 

The first step of the code(s) is the creation of region boundary nodes. Each segment’s 
length is computed, and nodes are then created along the region boundary based upon 
the user input nodal spacing. It is important to note that region boundaries should 
be created in a manner such that all common boundaries have the same vertices as 
endpoints. This boundary creation method insures that common boundary nodes will 
be created at the same locations between vertices. 
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1 2 3 4 5 

= Common nodes 


Figure 2.2: Creation of boundary nodes. 

Due to the method by which the region boundaries are defined, the nodal “com- 
pression” algorithm is straightforward. The code creates a list of compressed nodes, 
separate from the “uncompressed” boundary node list, beginning with the unique 
boundary nodes from the first region. The code then compares the boundary nodes 
of the remaining regions to the compressed node list, adding only uniquely defined 
nodes to the compressed list (see Figures 2.2 and 2.3). Additionally, the compression 
algorithm tracks and maintains a list of global node information: PEC TM, PEC TE, 
and RES global node numbers and their corresponding values. No nodal compression 
is necessary if the geometry contains a single region. 

Nodal Creation Interior to Domain Boundaries 

The method used for internal node creation can be summarized by the following steps 
[7] (see Figure 2.4): 
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1. Determine the minimum and maximum Y value of the region. 

2. Adjust the average element spacing for the region based upon the Y m j n 
and 7 ma x values of the region and the user-input nodal spacing. 

3. Make imaginary horizontal cuts at Y values determined from 1) and 2) 
above. 

4. For each horizontal cut, determine the intersection point X coordinates 
by cycling through region boundary segments and determining which seg- 
ments are intersected. Each boundary segment has endpoints (XI, Yl) 
and (X2,Y2). The imaginary line drawn at Y=H intersects a straight line 
segment if: 

(a) (Yl - H )(Y2 - H) < 0 or 

(b) (Yl - H)(Y2 -H) = 0 and (H > Yl or H > Y2) 

If intersection of the imaginary horizontal line with a straight line segment 
occurs, the X-value of the intersection point can be found from: 


(h-y 1 )(x 2 -x 1 ) 

^ inter ^ + (y 2 _ y a ) 


(2.1) 


To determine whether or not intersection occurs between the imaginary 
horizontal line and an arc segment, consider the equation of a circle: 


(X - X c f + (Y - Y c f = R 2 (2.2) 

The intersection points of the imaginary horizontal line (Y=H) with an 
arc segment are then given by: 

Winter = ± s[& ~ ( H ~ Y? (2 3) 

Intersection occurs if: 

(a) R 2 -(H- Y c ) 2 > 0 and 

(b) The arc segment passes through Xj n ^ er given above. 
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5. Node creation is considered on the horizontal line of 3) between intersec- 
tion points found in 4). Since the code that uses Delaunay triangulation 
is designed to mesh only convex polygonal regions, just two intersection 
points will be found in 4). The code that uses the advancing front method 
computes and considers 2n intersection points two at a time, with poten- 
tial nodes located between each set of points (see Figure 2.4). Whether or 
not a node is finally created is determined by whether or not the potential 
internal node is sufficiently spaced apart from all region boundary nodes. 
If (X,Y) is a potential nodal location, (Pi,Qi) is a boundary node, and 
C is a constant which is dependent on the average element size D of the 
region, usually defined as .ID, then a new node will be created if: 

(■ X - Pi ) 2 + (Y- Qi ) 2 < C 2 (2.4) 

The logic used to compute intersection points in 4) above has the disadvantage 
that intersection points will be computed for each segment, possibly resulting in a 
redundant list of points. Thus, some “post processing” is implemented to sort out 
the correct pair(s) of intersection points between which new nodes may be created. 

When nodal creation is complete for a region, the list of X and Y coordinates for 
nodes in the region are passed to the triangulation section of the code. 

2.4.4 Triangulation 

Delaunay Triangulation 

This section presents Watson’s [5] algorithm as described and implemented in [6]: 

1. Sort the N points to be triangulated in ascending sequence of their x coordinate. 

2. Define the vertices of the “supertriangle” in a counter-clockwise manner. The 
supertriangle is the initial Delaunay triangulation consisting of a single triangle. 
The size and shape of the supertriangle can be arbitrary, but should contain 
all of the points to be triangulated. The supertriangle should begin a list of 
triangles formed and flagged as incomplete. 
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3. Introduce a new point (X new , Y new ) from the list of sorted points. 

4. For each incomplete triangle, perform steps 5 to 9. (see Figure 2.6 for an 
illustration of terms defined in steps 5-9.) 

5. Compute the coordinates of the triangle’s circumcentre (X Ci Y c ) and the square 
of the triangle’s circumcircle radius R 2 . An illustration of a triangle’s circum- 
circle radius is provided in Figure 2.5. 

6. Compute the square of the distance in the x direction from the new point to 
the triangle’s circumcentre: D x 2 = ( X c — X ncw ) 2 . 

7. UD X 2 > R 2 then the circumcircle of the triangle cannot be intersected by any 
of the remaining points. Flag this triangle as complete and skip steps 8 and 9. 

8. Compute the square of the distance from the new point to the triangle’s cir- 
cumcentre: D 2 = D x 2 + (y c — Vnew ) 2 • 

9. If D 2 < R 2 then the new point intersects the circumcircle of the triangle. Re- 
move this triangle from the list of triangles formed and store three sets of vertices 
defining the triangle’s edges on a list of edges. If D 2 > R 2 > the new point lies 
outside the triangle’s circumcircle and the triangle remains unmodified. 

10. Search the list of edges and remove all edges which occur twice in the list. 
This step will remove all edges which are interior to the polygon formed by the 
intersected triangles. 

11. Form new triangles with vertices in a counter-clockwise fashion using the new 
point combined with the vertices of each edge of the list of edges (see Figure 
2.7). Flag each new triangle as incomplete. 

12. Repeat steps 3-11 until all remaining nodal points have been introduced. 

13. Complete the Delaunay triangulation by removing all triangles which contain 
one or more of the supertriangle vertices. 
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Figure 2.5: Illustration of circumcircle radius. 



• = (Xnew,Ynew) 

□ = triangle circumcentre (Xc,Yc) 

o = nodal points 
R= circumcircle radius 

Figure 2.6: Illustration of Delaunay triangulation geometry. 
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List of edges: 
AB,BC,CA 

Incomplete triangles: 
abc,abd,bcd,cad 



Figure 2.7: Delaunay triangulation: Formation of new triangles. 



initial generation front Cl = best point C best point C chosen 

AB - current segment AB de-activated AB de-activated 

suitable points Cl & C2 BC de-activated AC activated 

AC activated CB activated 

■ = current active segment 

Q = generation front 

Figure 2.8: Generation front illustration. 


Advancing Front Triangulation Implementation 

This section will describe details needed for numerical implementation of an advancing 
wavefront triangulation algorithm based upon a classic paper by Lo [7]. 

A term that needs introduction is “generation front”. In the advancing wavefront 
method, the generation front is the collection of line segments which are available to 
form a new triangular element. The generation front changes continuously as each 
new element is created, (see Figure 2.8). Triangulation terminates when no segments 
remain in the generation front. 

Steps that need to be executed prior to triangulation are described by the follow- 
ing: 
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1. Define the generation front for the current region. Initially it is the col- 
lection of boundary segments of the region listed in a counter clockwise 
fashion (see Figure 2.8). 

2. Activate the generation front segments. Activation of the segments amounts 
to defining a flag to indicate which segments of a list of segments are part 
of the generation front. 

3. Define the list of nodes to be triangulated. Maintain a link between local 
(referenced to the current region) node numbers and global (referenced 
to the entire geometry) node numbers. Nodes are classified as generation 
front nodes or interior nodes. A point that is suitable for triangulation 
with the current segment must be either a generation front or an interior 
node, eliminating the possibility of overlapping elements. 

Triangulation then proceeds as follows: 

1. Define the current active segment (CAS), segment A-B, from the list of 
generation front segments. 

2. Find a set of potential triangulation points that have the smallest norms. 
The size of this set of points is dependent upon the user and is set as a 
program parameter. If C is a potential triangulation point, then the norm 
is written as [AC 2 -f CB 2 ], and defined by: 

[AC 2 +CB 2 ] = [(. X A -XcfHyA-ycf] + [{X B -XcfHyB-Yc ) 2 ] (2.5) 

3. Since the generation front is defined in a counterclockwise manner, all valid 
triangulation points will be situated to the left of the CAS. Points found 
from the previous step are flagged as valid or invalid based upon whether 
they are to the left or to the right (or co- linear) to the CAS, respectively. 
The two points to the left of the CAS with the smallest norms, defined as 
C\ and C 2 , are considered the best potential triangulation points and are 
compared in the next step. 
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a=. 1443 a=.125 a=.1102 a=.1083 

Figure 2.9: Triangle quality examples. 

4. The quality of triangles C\ and C 2 are computed from the following rela- 
tions, where A ABC = the area of triangle ABC : 


AABCi 

Ql “ AB 2 + BC\ + Ci A 2 

(2.6) 

A C x BC 2 

CiB 2 + BC\ + CiC\ 

(2.7) 

A^lCiC^ 

1 _ AC\ + CiC\ + C 2 A 2 

(2.8) 

A a = max(/3i, 8 1 ) 

(2.9) 

aabc 2 

02 “ AB 2 + BCi + C 2 A 2 

(2.10) 

A C 2 BCi _ 

C 2 B 2 + BCI + CiC\ Pl 

(2.11) 

&AC 2 C\ £ 

~ AC\ + C 2 C\ + Cl A 2 ~ 1 

(2,12) 

A 2 = mai(/3 2 , 8 2 ) 

(2.13) 


Ci is chosen if > a 2 A 2 , otherwise C 2 is chosen. Some examples of 

computed triangle quality are shown in Figure 2.9. The larger the a value, 
the more suitable the element for the finite element method. Details of 
triangle quality criterion given by the above equations should be referenced 

to [7]. 
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5. To insure that the new element ABC does not overlap an existing ele- 
ment, it should be verified that the most suitable point C is either on the 
generation front or an interior point. 

6. Add the new element to the list of elements for the current region. 

7. De-activate the current edge. 

8. Activate and de-activate segments accordingly. The three cases that can 
occur are illustrated in Figures 2.8 and 2.10 and are described as follows: 

(a) De-activate old segments AB and BC, activate new segment AC. 

(b) De-activate old segment AB, activate new segments AC and CB. 

(c) De-activate old segments AB, BC, and CA (corner element). 

The logic used for the appropriate activation and de-activation of segments 

is illustrated with the following pseudo-code: 

De-activate current segment AB 
If segment AC is already active, flag AC as old 
If segment CB is already active, flag CB as old 
If AC is not old then 
Activate AC 

Add AC to active segment list 
end 

If CB is not old then 
Activate CB 

Add CB to active segment list 
end 

If one or more new segments are added then 
Point to last added segment 
If only one new segment is added then 
De-activate old segment 
Decrement active segment counter 
end 

else (corner element logic) 

Find new active segment by decrementing generation 
front pointer 
end 


9. Find the new CAS by pointing to the last active segment in the generation 
front segment list. 
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corner element 
encountered 


AC de-activated 
BC de-activated 
CA de-activated 
new active segment found 


■ = current active segment 


= generation front 


Figure 2.10: Corner element logic. 


10. If three or more active generation front segments remain, return to step 2) 
and continue the triangulation process. 


2.5 Meshed Geometries 

Three meshed geometries are shown in Figures 2.11 through 2.13. The meshes cor- 
respond to a PEC cylinder, two PEC cylinders, and a semicircular PEC cylinder 
with dielectric skin. The files used to create the geometries are given in Figures 2.14 
through 2.16. 
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Figure 2.11: Meshed single cylinder geometry. 
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Figure 2.12: Meshed dual cylinder geometry. 
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Figure 2.13: Meshed semicircular cylinder with material skin geometry. 
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5 

0 70 
0 60 
0 0 
0 -60 
0 -70 
2 
4 

’line’ 2 1 'none’ 
'arc’ 1 5 ’none* 

3 'ccw' 

'line’ 5 4 ’none’ 
'arc' 4 2 ’ pectin ’ 

3 ’ cw’ 

4 

'line' 1 2 ’none’ 
'arc' 2 4 ’pectin' 
3 ’ cw’ 

’line’ 4 5 ’none’ 
’arc’ 5 1 ’none’ 

3 ’ccw’ 


Figure 2.14: Single PEC cylinder geometry input file. 
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10 

-11 22 
-10 22 
0 22 
10 22 
11 22 
-11 -22 
-10 -22 
0 -22 
10 -22 
11 -22 
3 
8 

’line’ 10 5 ’none' 
’line’ 5 4 ’none’ 
’arc’ 4 2 ’ pectin ’ 


3 ’cw’ 



’line’ 2 

1 

'none ' 

’line’ 1 

6 

'none ' 

’line’ 6 

7 

'none' 

’arc’ 7 9 

> 

pectm' 


8 ’cw’ 

’line’ 9 10 ’none’ 
4 

’arc’ 5 1 ’none’ 

3 ’ccw’ 

’line’ 1 2 ’none’ 
’arc’ 2 4 ’pectm’ 

3 ’cw’ 

’line’ 4 5 ’none’ 

4 

’line’ 10 9 ’none’ 
’arc’ 9 7 ’pectm’ 

8 ’cw’ 

’line’ 7 6 ’none’ 
’arc’ 6 10 ’none’ 

8 ’ccw’ 


Figure 2.15: Dual PEC cylinder geometry input file. 
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11 

0.000000 0.000000 
0.000000 0.500000 
0.000000 0.829775 
0.000000 6.850000 
0.000000 12.882500 
0.000000 13.200000 
0.000000 13.700000 
42.931801 12.882500 
45.237396 13.200000 
45.737396 12.840371 
45.737396 13.700000 
4 
7 


’line’ 

10 

» 11 'none 

’line ’ 

11 

7 

'none ' 

’line ’ 

7 

6 

'none' 

’line’ 

6 

9 

'none' 

’line’ 

9 

2 

'none' 

’line’ 

2 

1 

'none' 

’line’ 

6 

1 

10 

' none ' 

’line’ 

2 

9 

'none' 

’line’ 

9 

6 

'none' 

’line’ 

6 

5 

'pectin' 

’line’ 

5 

8 

'none' 

’line’ 

8 

3 

'none' 

’line’ 

3 

3 

2 

'pectin' 

’line’ 

8 

5 

' none ' 

’line ’ 

5 

3 

'pectin' 

’line’ 

4 

3 

8 

'none' 

’line’ 

6 

7 

'none' 

’arc’ 

4 ’ccw 

7 1 

) 

none' 

’line’ 

1 

2 

' none ' 


'arc* 2 6 ’pectin* 

4 ’cw’ 

Figure 2.16: Semicircular PEC cylinder geometry with dielectric skin geometry file. 
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Chapter 3 

Genetic Algorithms 


Genetic algorithms (GAs) [3], as applied to numerical optimization problems, attempt 
to imitate the concept of “survival of the fittest” in nature by simulating the processes 
of genetic reproduction, recombination, mutation, and evolution. A description of the 
basic concepts and terms related to GAs as well as some simple applications will be 
presented in the following sections. 

3.1 Basic Theory and Terms 

Each parameter of a problem of interest is constrained to fall within a specified range. 
The range for each parameter is discretized into a set of equally spaced values. The set 
size is defined using a power of two, allowing each parameter value to be encoded as a 
binary sequence. Each of these binary sequences is termed a gene. As an illustration 
of a gene, see Figure 3.1. In this example, the parameter that is discretized is the 
separation distance between a PEC plate and a resistive card. The upper and lower 
bounds of this parameter are defined as 8 cm and 1 cm, respectively, with the chosen 
discretization being defined by 2 3 or 8 values. Thus, the granularity for this parameter 
is 1 cm. In theory, a GA parameter can be defined using any range, and discretized 
using any power of two. 

Each of n genes of a set are appended together to form a chromosome. As an 
illustration, consider the geometry shown in Figure 3.2. A chromosome is simply a 
single binary number encoded to represent a specific value of each of the n parameters 
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PEC 


R-card 


d 

* ► 


--H-H-H-j 

1 cm 8 cm 


parameter: separation distance d 
min: 1 cm 
max: 8 cm 

discretization: 8 (2 A 3) 


gene code 

parameter 

000 

1 cm 

001 

2 cm 

010 

3 cm 

011 

4 cm 

100 

5 cm 

101 

6 cm 

110 

7 cm 

111 

8 cm 


Figure 3.1: Illustration of a gene. 


simultaneously. In this example, the two parameters represented by a chromosome 
are separation distance and resistive sheet value. 

Once the concepts of gene and chromosome have been introduced, it is useful to 
reference a flowchart similar to that of Figure 3.3 to help describe the main processes 
of a GA. A collection of chromosomes is termed a population. Initially, a population is 
created by generating a random collection of chromosomes. Once an initial population 
is created, the “fitness” of each member chromosome is computed. A chromosome’s 
fitness, in general, can be thought of as a numerical assessment of some user-defined 
“quality” measure based upon the given set of input parameters. After an initial 
random population has been generated, the creation of a new population begins. 
The original population is altered by implementing, in a probabilistic manner, the 
concepts of crossover, mutation, and reproduction. 

Crossover involves selecting, based on fitness, two chromosomes or “parents” of the 
current population and combining some part of the first parent’s string information 
with the complementary part of second parent’s string information. A crossover point 
is randomly selected, with a pair of offspring resulting from “splicing” the binary 
information from both parents together. See Figure 3.4. 
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PEC 


R-card 


parameter 1 : seperation distance d 
min: 1 cm 
max: 8 cm 

discretization: 8 (2 A 3) 


d 


parameter 2: R-card value RES 
min: 100 QAd 
max: 800 £2/d 
discretization: 8 (2 A 3) 


gene 1 

d value 

gene 2 

RES value 

000 

1 cm 

000 

100 O/D 

001 

2 

001 

200 

010 

3 

010 

300 

011 

4 

on 

400 

100 

5 

100 

500 

101 

6 

101 

600 

110 

7 

110 

700 

111 

8 

111 

800 


Chromosome examples: 

010 +110 = 0101 10 : d=3, RES=700 
101 + 001 = 101001 : d=6, RES=200 
111+000= 1 1 1000 : d=8, RES=100 

\/ t 

genes chromosomes 


Figure 3.2: Illustration of a chromosome. 


Mutation involves the toggling of a small percentage of bits in the population 
of chromosomes. Increasing the number of mutations allows a GA more freedom to 
search for solutions outside the present population. Goldberg [9] states: 


That the mutation operator play a secondary role in the simple GA, we 
simply note that the frequency of mutation to obtain good results in 
empirical genetic algorithm studies is on the order of one mutation per 
thousand bit (population) transfers. Mutation rates are similarly small 
(or smaller) in natural populations, leading us to conclude that mutation 
is appropriately considered as a secondary mechanism of genetic algorithm 
adaptation... when used sparingly with reproduction and crossover, it is 
an insurance policy against premature loss of important notions. 

Reproduction involves selecting an individual chromosome based upon its fitness 
as compared to others within the current population. Once this chromosome has been 
selected, it is merely copied and placed into the new population. 

Crossover, mutation, and reproduction continue until a new population is com- 
pleted. A complete population is termed a generation. After each new generation 
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Figure 3.3: Flowchart of a conventional genetic algorithm. 
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parent #1: 10101 0010101 
parent #2: 11100 1100101 

offspring #1 : 10101 1100101 
offspring #2: 11100 0010101 

crossover point 

Figure 3.4: Illustration of crossover. 

is created, it is analyzed to determine if a specified termination criterion has been 
satisfied. If the criterion have been satisfied, the algorithm is complete, otherwise 
new generations continue to be created. 

3.2 Differences between Genetic Algorithms and 
Traditional Methods 

Goldberg cites four main ways in which GAs are different from normal optimization 
and search procedures: 

1. GAs work with a coding of the parameter set, not the parameters them- 
selves. 

2. GAs search from a population of points, not a single point. 

3. GAs use payoff (objective function) information, not derivatives or other 
auxiliary knowledge. 

4. GAs use probabilistic transition rules, not deterministic rules. 
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The advantage of working with a coding of the parameter set is described using 
the concept of the “black box”. Encoding the parameter set to form a single binary 
number is much like using a set of binary switches as input to a black box which 
has a single output. The input switches are toggled to “on” or “off” in an organized 
fashion (initially at random), producing an output value representing the fitness of the 
particular binary string. This method of using only the single black box output value 
allows GAs to be removed from the problems associated with the internal workings 
of the box. (Some examples of the “problems” that may be associated with the black 
box are continuity, existence of derivatives, etc.) 

An important example that illustrates the advantages of using information from a 
population of points, as opposed to information from a single point, is a search space 
that contains many peaks. This multi-peak situation is where many optimization 
methods will locate false peaks, due to their lack of global information. Since GAs 
use information from many points simultaneously, their probability of finding an 
incorrect peak is greatly reduced in comparison with point to point methods. 

Another characteristic of GAs that makes them desirable over traditional methods 
is the fact that are essentially independent of auxiliary information. They rely only 
upon the previously discussed fitness information. Goldberg points out that: 

... the refusal to use specific knowledge when it does exist can place an 
upper bound on the performance of an algorithm when it goes head to 
head with methods designed for that problem. 

One final difference discussed here is that GAs use probabilistic methods. GAs 
initially use random search to create a broad collection of candidate solutions. Once a 
variety of candidate solutions exist, GAs use crossover, reproduction, and mutation to 
lead the search towards unchartered regions of the solution space where improvement 
is likely. This form of search guidance is possible due to the fact that each new 
population helps to bias the search based upon what has worked best so far. A 
search that guides itself based upon known information is more effective than a purely 
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random search which may contain a very large number of possible solutions, and has 
no specific direction to its search. 

3.3 Example Applications 

3.3.1 Infinite Geometry 

As an example of an application of a GA, consider a TM polarized plane wave incident 
upon an infinite PEC plate with resistive sheet geometry as shown in Figure 3.5. 
Suppose a search is desired to find the best combination of resistance and separation 
distance to minimize the reflected field. Keeping this goal in mind, the experiment 
proceeds as follows. Let the frequency of interest be 0.3 GHz such that the free 
space wavelength is 1 m. Additionally, let d be the separation distance between 
the PEC and the resistive sheet in m, and RES be the resistive value of the sheet 
in Q/U. Initially the minimum, maximum, and granularity of d are set to .1, 1.0, 
and 2 6 — 64 parameter levels (Ad ~ 1.43cm), respectively. Likewise, the minimum, 
maximum, and granularity of RES are set to 10, 1280, and 2 7 = 128 parameter 
levels (A RES 10f2/D), respectively. The total number of possible parameter 
combinations in this example is thus 2 13 = 8192. 

A search was performed, using the parameters described above, using the GA 
software GENESIS [4]. (Specific details about the use of GENESIS can be found in 
[10].) In general, a setup program is provided that acquires information about each 
search. The input information is recorded by two files, “in.esl” and “template. esl”, 
where the extension “esl” can be any user defined extension name. The template file 
contains the pertinent information about the run, and is described below. 

Total Trials Controls the number of times a numerical assessment of a chromosome 
is calculated. 

Population Size Specifies the number of chromosomes per population. 

Structure Length Specifies the number of bits per chromosome. 
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Infinite PEC plate 



Infinite resistive sheet, 
value=RES 


Figure 3.5: Geometry corresponding to a two parameter search for optimum distance 
of separation and resistive sheet value. 

Crossover Rate Specifies the percentage of each population (number of genes in 
each population) in which crossover occurs. 

Mutation Rate Specifies the percentage of bits which are mutated (toggled). 

Generation Gap Specifies the fraction of the population which is replaced in each 
generation. 

Scaling Window A brief explanation of the scaling window is given as follows: 

The fitness u(x) of a structure x can be defined by: 

u(x) = f max - f(x) (3.1) 

where f max is the maximum value that f(x) can take on in the given search 
space. Setting the scaling window W to a value greater than 0 sets / max to the 
greatest value of /(x) that has occurred in the last W generations. 
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Experiments 
Total Trials 
Population Size 
Structure Length 
Crossover Rate 
Mutation Rate 
Generation Gap 
Scaling Window 
Report Interval 
Structures Saved 
Max Gens w/o Eval 
Dump Interval 
Dumps Saved 
Options 
Random Seed 
Rank Min 


1 

1000 

50 

13 

.6 

.001 

1 

5 

100 

10 

2 

0 

0 

bCefgo 

123456789 

.75 


Figure 3.6: “in” file contents corresponding to a coarse two parameter search for 
optimum separation distance and resistive sheet values. 


Report Interval Specifies the number of trials that occur between collection of ex- 
periment information data. 

Structures Saved Specifies the number of best structures that are to be saved and 
placed in the “.min” file which is described later. 

Max Gens w/o Eval Sets the number of consecutive generations which are per- 
mitted without any fitness evaluation. 

Options Control the output and program execution which are documented in [10]. 

Rank Min Specifies the minimum number of offspring used in ranking and is used 
only if the appropriate option is set. 


The contents of the “in” and “template” files for the case described above are 
shown in Figures 3.6 and 3.7. The user creates a subroutine to assign each param- 
eter combination a fitness value. In this example, the magnitude of the reflected 
field represents the fitness of the particular combination of separation distance and 
resistance. 
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genes: 2 


gene 0 
min: 0.1 
max : 1 
values: 64 
format: '/,10.4f 

gene 1 
min: 10 
max: 1280 
values: 128 
format: */,10.4f 

Figure 3.7: “template” file contents corresponding to a coarse two parameter search 
for optimum separation distance and resistive sheet values. 

Once execution of the search has terminated, a “min” file is created. This file 
contains the best parameter combinations encountered, each combination’s fitness, 
and the generation and trial counters at the time of the first occurrence of the com- 
bination. The results from this run are summarized in Table 3.1. 

Basic transmission line theory helps to explain the results of this search, as shown 
in Figure 3.8. Using impedance matching concepts, a short circuit is converted to 
an open circuit when transformed using a transmission line whose length is an odd 
multiple of A 0 /4 at the design frequency / 0 . Once the short circuit is converted to an 
open circuit, setting R a equivalent to the characteristic impedance of the line (in this 
case rjo = y^o/eo ~ 377Q) is sufficient to eliminate reflection. Note that, in contrast 
to point to point methods, two solutions were found in this search space. 

An important observation should be made at this point. When discretizing con- 
tinuous input parameters, it is possible that the discretization will not include a 
parameter level that could potentially have led to a combination with better fitness. 
For example, the fitness values of the two solutions found in this search are very 
similar, but suppose that levels of 30 dB down are not sufficient. Lower levels may be 
possible, but may not have been found due to better parameter values being missed 
by the search even if they fall between the max and min of the defined parameter 
window. As another experiment, it is useful to narrow the parameter search windows 
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Table 3.1: Results of coarse two parameter search for optimum distance of separation 
and resistive sheet value. 


d 

(m) 

R. 

(ft/D) 

Magnitude 

(dB) 

generation 

trial 

0.7429 

380.0000 

-3.2809e+01 

11 

278 

0.2571 

370.0000 

-3.2348e+01 

14 

317 

|[ 


-3.2809e+01 

7 

201 

0.2429 

370.0000 

-3.2348e+01 

25 

403 

0.7571 

380.0000 

-3.2809e+01 

9 

236 

0.7571 

390.0000 

-3.0924e+01 

8 

212 

0.7429 

390.0000 

-3.0924e+01 

9 

241 

0.2571 

390.0000 


2 

Oil 

0.2429 

390.0000 


6 

179 

0.2571 

380.0000 

-3.2809e+01 

4 

123 
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Figure 3.8: Basic transmission line theory. 
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genes : 2 


gene 0 
min: 0.1 
max : 0.5 
values : 64 
format: '/,10.4f 

gene 1 
min: 10 
max : 640 
values: 128 
format: */,10.4f 

Figure 3.9: “template” file corresponding to a finely discretized two parameter search 
for optimum separation distance and resistive sheet value. 

while keeping the number of discretizations constant. This narrowing allows finer 
parameter values to be tested based upon combinations that were promising in the 
previous search. 

As an illustration, the parameter window for d is modified to have minimum and 
maximum values of 0.1 and 0.5, respectively, with granularity of 2 6 = 64 parameter 
levels (Ad ~ ,635cm). The parameter window for RES is also narrowed to have 
minimum and maximum values of 10 and 640, respectively, with granularity of 2 7 = 
128 parameter levels (A RES ~ 4.96U/D). The narrowing of the parameter windows 
in this case allows further investigation of the solution near 0.25 m and 380 fi/D. 
The “template” file for this example is shown in Figure 3.9. The “in” file remains 
the same as in the first example. The results of this search are shown in Table 3.2. 
Utilizing a more narrowly defined search space, the GA is able to focus upon the 
minimum that occurs at d = A 0 /4 and RES = 770 ~ 377ft. 


3.3.2 Semi-infinite Geometry 

The previous calculations were performed using structures of infinite extent. A more 
practical geometry to apply a GA to is that of Figure 3.10. A two dimensional FEM 
analysis is used to calculate the backscatter levels of various separation and resistive 
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Table 3.2: Results of a finely discretized two parameter search for optimum distance 
of separation and resistive sheet value. 



R, 

(f2/D) 

Magnitude 

m 

generation 

trial 



-4.0107e+01 

7 

175 

0.2460 

372.1260 

-3.7083e+01 

4 

120 

0.2524 

382.0472 

-3.9960e+01 

3 

93 



-3.8083e+01 

10 

242 



-3.6953e+01 

7 

178 



-4.2519e+01 

5 

143 

0.2460 

367.1654 

-3.4868e+01 

7 

168 

11 


-3.6399e+01 

5 

146 

0.2460 

387.0079 

-3.4800e+01 

2 

74 

0.2524 

387.0079 

-3.6401e+01 

14 

274 
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Figure 3.10: Second geometry corresponding to a two parameter search for optimum 
distance and resistive sheet value. 


sheet values with structured data scan results shown in Figure 3.11. A two dimen- 
sional GA search was performed using the parameter space defined by the “template” 
file contents shown in Figure 3.12. As expected, the peak was found with optimum 
parameter combinations given by the “min” file shown in Table 3.3. The peak levels 
found in the search space are not as low as obtained through the use of the structured 
grid. This can be contributed to discretization of continuous parameters. The peak 
computed by the scanned parameter set resulted from a parameter combination that 
was not possible based upon the parameter discretizations that were defined for this 
search. 


3.3.3 Conclusions 

The previous examples represent searches for optimum combinations of just two pa- 
rameters. In these cases, it is more reasonable to search the solution space and 
describe the data using an illustration similar to that of Figure 3.11. There may be 
times where this is not practical for two reasons. 
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Figure 3.11: 2-D scan for minimal reflected field from a 5 inch PEC plate with one 5 
inch resistive card. 0° incidence, TM polarization, / = 1 GHz. 


genes: 2 

gene 0 
min: 100 
max: 400 
values: 32 
format: '/,10.4f 

gene 1 
min: 0.1 
max : 0.5 
values : 32 
format: '/,10.4f 

Figure 3.12: “template” file corresponding to the search for optimum resistive sheet 
and separation distance values of the finite geometry. 
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Table 3.3: Search results for optimum resistive sheet and separation distance values 
of a finite geometry. 


R. 

(«/□) 

d 

w 

Echo Width 
(dBm) 

generation 

trial 

225.8065 

0.1903 

-3.3205e+01 

7 

180 

196.7742 

0.1903 

-2.9238e-f 01 

7 

179 

225.8065 

0.2032 

-3.4312e+01 

4 

134 

235.4839 

0.2032 

-3.1559e+01 

4 

131 

235.4839 

0.1903 

-3.0283e+01 

0 

42 

254.8387 

0.2032 

-2.6644e+01 

3 

103 

216.1290 

0.2032 

-3.4452e+01 

6 

159 

245.1613 

0.1903 

-2.7747e+01 

3 

94 

216.1290 

0.1903 

-3.4690e+01 

5 

141 

245.1613 

0.2032 

-2.8836e+01 

6 

165 


44 




The first reason is that the computational time to create a set of data may be 
quite large. For example, if it is desired to create a data set of 100 distances and 
100 resistance values, this represents 10,000 data points. If each data point takes 1 
second to compute, this amounts to close to three hours of computation time. In this 
case it would be much more desirable to use a GA to search the solution space using 
a fraction of the number of data points. 

A second and more important reason that it may not be practical to perform a 
parameter scan would be when more than two parameters are of interest in a search. 
GAs are designed in general such that they search the solution space of a single 
parameter problem using the same algorithm that they search the solution space of 
a multi-parameter problem. 
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Chapter 4 

Scattering Reduction Techniques 


The basic components of the geometry of interest are illustrated in Figure 4.1. The 
geometry consists of a conducting bulkhead, a dielectric skin, and a resulting internal 
region. The goal of this chapter is to provide an investigation of the scattering behav- 
ior of each component of the geometry separately before attempting to understand 
the interaction of all the components simultaneously. 

4.1 Conducting Bulkhead Characterization 

The conducting bulkhead of the scatterer of interest is modeled using various geome- 
tries shown in Figure 4.2. It is desired to reduce the backscatter levels of the PEC 
bulkhead(s) over a range of frequencies and incidence angles. The backscatter charac- 
teristics of each of these geometries are shown in Figure 4.3. Note that a conducting 
bulkhead with a wedge shaped PEC termination has lower backscatter levels over the 
frequency range of interest. 
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Figure 4.2: Bulkhead geometries. 
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Figure 4.3: Echo widths for the planar, circular, and wedge shaped bulkheads. 
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Figure 4.4: Infinite PEC wedge geometry. 

4.2 Conducting Bulkhead Shape Considerations 

4.2.1 Wedge Angle 

In the previous section, it is illustrated that a wedge shaped conducting bulkhead with 
a wedge angle of approximately 15° has lower backscatter levels than both the circular 
and planar bulkheads. To further investigate the dependence of the diffracted field on 
wedge angle, consider the geometry shown in Figure 4.4. Using a GTD approximation 
[11], the diffraction characteristics of this geometry are computed and given in Figure 
4.5. Additionally, Figure 4.6 illustrates the frequency response of an infinite PEC 
wedge. Figure 4.5 suggests that TM polarization limits how much of an effect wedge 
angle can have on the diffracted field. The data also suggests that variation of wedge 
angle has little effect on the diffracted field when the wedge angle considered remains 
small. 


4.2.2 Planar Bulkhead Orientation 

Consider the geometry shown in Figure 4.7. Figure 4.8 illustrates how the echo width 
level changes as the angle of incidence is varied from 0° to 30°. Note the maximum 
echo width level is approximately 0 dBm. This level changes with frequency, as the 
electrical size of the bulkhead changes. The maximum level can be used as a reference 
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0 30 60 90 

Wedge Angle (deg) 

Figure 4.5: Magnitude of the diffracted field of an infinite PEC wedge, / = 3 GHz 
0° incidence. 
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Figure 4.6: Echo Widths for an infinite PEC wedge, 15° wedge angle, 0° incidence. 





Figure 4.7: Planar bulkhead geometry illustrating angle of incidence. 

when considering the effects of rotating a conducting planar bulkhead with dielectric 
skin. 

To study the tilt angle 0 T of the planar bulkhead with dielectric skin, consider 
Figure 4.9. The dielectric skin has an approximate thickness of .3175 cm with e Ti ,k tn = 
4. Figures 4.10 and 4.11 illustrate backscatter echo widths for rotation angles of 
0r=lO°, 20°, and 30° for TM and TE polarizations, respectively. The echo width 
levels are similar to the levels of the bulkhead in free space. The maximum echo 
width level of the rotated bulkhead geometry returns to the maximum level of the 
bulkhead in free space for rotation and incidence angles in the range 0° to 30°. This 
may indicate that rotation of the bulkhead is not an effective method of reducing 
backscatter. 

4.2.3 Bulkhead Junction Blending 

In the first section of this chapter, it was shown that the wedge shaped conducting 
bulkhead had lower echo width levels over the frequency range of interest than did 
the planar or circular conducting bulkheads. These lower levels are a result of the 
wedge shape bulkhead providing smoother shape discontinuities than the planar and 
circular geometries. 

Based upon these preliminary results, three junction blending techniques are con- 
sidered in shaping the junction between the wedge and the bulkhead as shown in 
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Ech 



Figure 4.8: Backscatter echo widths of the planar bulkhead with respect to angle of 
incidence, / = 3 GHz. 



Figure 4.9: Planar bulkhead geometry illustrating angle of rotation 6 T . 





Figure 4.10: Echo widths of rotated planar bulkhead geometry with TM polarization. 
/ = 3 GHz. 0 r = lO° - solid, 20° - long dashed, 30° - short dashed. 



Figure 4.11: Echo widths of rotated planar bulkhead geometry with TE polarization. 
/ — 3 GHz. 0 T — 1 0° - solid, 20° * long dashed, 30° - short dashed. 
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GEOMETRIES ENTIRELY PEC 


30.5 cm- 


45.2 cm- 



circularly blended termination 
linearly blended termination 
linear wedge termination 



circularly blended 
linearly blended 


Figure 4.12: Three blended termination types. 


Figure 4.12. The geometry names that are used to differentiate between the three 
terminations tips are the circularly blended termination, the linearly blended termi- 
nation, and the linear wedge termination. The geometries consist of a circular end 
cap, a rectangular body, and a termination tip. In each case, the body of the bulkhead 
is elongated so that scattering from the tip, junction, and endcap can be isolated and 
studied using time domain analysis. 

Swept frequency calculations for each of the three geometries are shown in Figures 
4.13 through 4.15. The three scattering mechanisms that appear as ripples in the 
transient scattering signature(s) are the front termination (~-3 ns), the bulkhead 
junction (0 ns), and the end termination (~2-4 ns). 

Figures 4.16 through 4.18 contain the extracted results of the junction alone. 
The extracted results are useful as they provide frequency domain information of 
the bulkhead junction independent of other scattering mechanisms. The linearly and 
circularly blended junctions are more effective in providing lower scattering levels over 
the frequency range of interest than the linear wedge junction. This is due to the 
circular and linear blended junctions providing much smoother transitions between 
the tip and junction than the linear wedge geometry. 
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Figure 4.14: Transient and frequency domain scattering signatures for the linearly 
blended termination geometry with TE illumination. 
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Figure 4.15: Transient and frequency domain scattering signatures for the circularly 
blended termination geometry with TE illumination. 
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Figure 4.16: Extracted tip scattering levels for the linear wedge termination junction 
with TE illumination. 
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Figure 4.17: Extracted tip scattering levels for the linearly blended termination junc 
tion with TE illumination. 









45.2 cm 



Figure 4.19: Planar conducting bulkhead geometry. 

4.3 Material Considerations 

The previous section explored methods of altering the conducting bulkhead geometries 
to lower the scattered field levels. This section will address the effects of dielectric 
materials and resistive cards when combined with a conducting bulkhead. 

An illustration of the geometry for analysis is shown in Figure 4.19. As indicated 
by the figure, there are several components of the geometry that effect its scattering 
signature: the dielectric skin; the internal material loading; and the placement of 
resistive cards and impedance surfaces. Each of these components will be discussed 
in the following sections. 

4.3.1 Dielectric Skin 

Wedge Angle 

iFrom Figure 4.5 it was concluded that wedge angle does not significantly effect the 
magnitude of the diffracted field from a conducting wedge when the wedge angle 
remains small. When considering a penetrable dielectric tip, however, wedge angle is 
of greater importance as will now be demonstrated. 

Figure 4.20 illustrates a model of the dielectric skin that will be incorporated 
with the bulkhead geometries. The size and shape of the skin may be constrained by 
other considerations. For now the length of the wedge formed by the skin and the 
thickness of the skin will be held constant to investigate the effects of varying the 
wedge angle. The backscatter characteristics for wedge angles of 10°, 20°, and 30° 
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are shown in Figures 4.21 and 4.22 for TM and TE polarizations, respectively. The 
results indicate that it is desirable to maintain a small wedge angle when considering 
a penetrable termination tip to keep backscatter levels as low as possible over the 
desired frequency range. 

Polarization and Angle of Incidence 

Polarization and angle of incidence also have an effect on the return from a dielectric 
slab. Consider the geometry shown in Figure 4.23. Figures 4.24 and 4.25 illustrate 
the transmission and reflection characteristics of a dielectric slab as a function of 
angle of incidence for TM and TE polarization, respectively. 

The transmission coefficient of the slab is of particular interest, as it is an in- 
dication of the amount of energy that penetrates the slab. The results show that, 
for both polarizations, as the angle of incidence approaches 90° (grazing), the slab 
rapidly becomes less transmissive. This is an important point. Consider Figure 4.26. 
The dielectric skin consists of two thin dielectric slabs combined into a wedge shape. 
When considering incidence angles near 0°, (see Figure 4.23) with respect to the 
upper slab approaches 90° , with 6 t for the lower slab having an approximate value of 
75°. 

Frequency and Relative Dielectric Constant 

Figures 4.27 through 4.30 show the effects of varying frequency and dielectric constant. 
For TM polarization, the slab becomes less transmissive as the electrical thickness 
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Figure 4.21: Extracted tip scattering levels for dielectric skin geometry, 
TM polarization. q= 10° - solid, 20° - dotted, 30° - dashed. 
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Figure 4.22: Extracted tip scattering levels for dielectric skin geometry, 
TE polarization. q= 10° - solid, 20° - dotted, 30° - dashed. 
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Figure 4.23: Dielectric slab geometry. 



Figure 4.24: Magnitude of reflection and transmission coefficients for a .3175 cm thick 
dielectric slab in free space. / = 3 GHz, e Tia i a b = 4, TM polarization. 
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Figure 4.25: Magnitude of reflection and transmission coefficients for a .3175 cm thick 
dielectric slab in free space. / = 3 GHz, e T>s i a b = 4, TE polarization. 


.45.2 cm 



Figure 4.26: Dielectric skin represented as two individual dielectric slabs. 
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Figure 4.27: Magnitude of reflection and transmission coefficients for a .3175 cm thick 
dielectric slab in free space. e TtS i ab = 4, TM polarization, 75° incidence. 

of the slab is increased. The slab also becomes less transmissive when the relative 
dielectric constant of the slab is increased. For TE polarization with an incidence 
angle of 75°, the slab remains highly transmissive in the given frequency and relative 
dielectric ranges. If the incidence angle of the slab approaches 90° (grazing), the 
magnitude of the (TE) transmission coefficient will decrease more dramatically as 
frequency is increased, as well as when the relative dielectric constant of the slab is 
increased. 

The next section considers a geometry consisting of a wedge shaped dielectric 
skin with interior material loading combined with the conducting planar bulkhead 
geometry (see Figure 4.19). It is found that backscatter is significantly reduced by 
loading the region between the conducting bulkhead and the dielectric skin with a 
lossy dielectric material. 

4.3.2 Internal Loading Parameters 

A set of e T values of interest for use in internally loading the wedge geometry of 
Figure 4.19 are shown in Table 4.1. The dielectric skin is assumed to take on a value 
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Frequency (GHz) 

Figure 4.28: Magnitude of reflection and transmission coefficients for a .3175 cm thick 
dielectric slab in free space. e Ti ,i a b = 4, TE polarization, 75° incidence. 



Real part of relative epsilon 


Figure 4.29: Magnitude of reflection and transmission coefficients for a .3175 cm thick 
dielectric slab in free space. / = 3 GHz, TM polarization, 75° incidence. 
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Figure 4.30: Magnitude of reflection and transmission coefficients for a .3175 cm thick 
dielectric slab in free space. / = 3 GHz, TE polarization, 75° incidence. 

of er.jkm = 4 or 6. Each internal loading parameter was evaluated for both TM and 
TE polarizations with each of the two values of skin dielectric constant. It was found 
for TM polarization that a skin dielectric constant of e r ,»ifctn = 6 was more effective 
in lowering the backscatter levels over the frequency range of interest. The internal 
loadings with the lowest backscatter levels are illustrated in Figures 4.31 and 4.32. 

To further reduce backscatter levels at low frequencies, a 500ff impedance is ap- 
plied to the flat portion of the conducting bulkhead, with results shown in Figures 4.33 
and 4.34. Implementation of this impedance surface, however, may not be realistic. 

Originally, it was thought that internally loading the interior wedge region with 
the material with the largest loss tangent would be the most effective in reducing 
backscatter. The previous results do not indicate that this is the case. 

Two possible reasons why the material with the largest loss tangent was not the 
most effective material in reducing backscatter are the following. First, upon inspec- 
tion of Table 4.1, materials with the largest loss tangent also have the largest value 
of |e r |, which effects the transmission and reflection characteristics of the material. 
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Table 4.1: Estimated lossy foam parameters. 


Re(e r ) 

Im(e r ) 

Loss tangent 

23.3 

-16.3 

.70 

18.8 

-9.01 

.48 

14.9 

-5.24 

.35 

11.8 

-3.16 

.27 

7.50 

-1.20 

.16 

6.00 

-0.70 

.12 

4.80 

-0.40 

.08 



Figure 4.31: Backscatter characteristics of the simple wedge geometry with internal 
loading of e r = 14.9-j5.24, e riS ^ n = 6, 0° incidence. 
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Figure 4.32: Backscatter characteristics of the simple wedge geometry with internal 
loading of e r =11.8-j3.16, e T ,skin = 6, 0° incidence. 



Figure 4.33: TM Backscatter characteristics of the simple wedge geometry with in- 
ternal loading of e r =11.8-j3.16, e r ,»kin — 6, 0° incidence, with 500f) impedance on flat 
portion of conducting bulkhead. 
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Figure 4.34: TE Backscatter characteristics of the simple wedge geometry with inter- 
nal loading of e r =11.8-j3.16, e Tit kin = 6, 0° incidence, with 500J2 impedance on flat 
portion of conducting bulkhead. 

Secondly, the change in scattered field behavior caused by the loaded interior region 
is considered in the presence of the dielectric skin. The interior region’s effects on the 
scattered field in the presence of the dielectric skin may be not be the same as when 
the interior region is considered independently. 

4.3.3 Resistive Cards and Impedance Surfaces 

Another method that can be used to influence the scattered field levels of a structure 
like that of Figure 4.19 is the placement of complex impedance surfaces. The simplest 
impedance surface would be a resistive (purely real) surface. Two possible locations 
for placement of such a surface are the conducting bulkhead or the exterior of the 
dielectric skin. The next sections address the effects of utilizing these surfaces. 

Resistive Cards on External Surface of Dielectric Skin 

Figures 4.35 and 4.36 compare the backscatter characteristics of a constant resistive 
card on the surface of the dielectric skin with that of a tapered resistive card for TM 
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Figure 4.35: Backscatter characteristics of the simple wedge geometry, 0° incidence, 
TM polarization. 

and TE polarizations, respectively. The tapered resistive card is modeled to have a 
low fi/D value at the conducting portion of the structure, tapered to a value of a few 
thousand fl/d at the tip. 

A resistive card of constant value yields better performance than a tapered resistive 
card since the card of constant value allows less energy to penetrate the interior 
region. In contrast, a tapered resistive card provides less scattering at the wedge tip, 
but allows for greater internal scattering. 

To interpret why a resistive card placed on the external surface of the dielectric 
skin is less effective when considering TE polarization, an infinite resistive card model 
can be used. 

Infinite Resistive Card 

Figures 4.37 and 4.38 illustrate the reflection and transmission properties, respec- 
tively, for an infinite R-card with 90° (grazing at 0°) being perpendicular to the plane 
of the card. 
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Figure 4.36: Backscatter characteristics of the simple wedge geometry, 0° incidence, 
TE polarization. 
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Figure 4.37: Magnitude of the reflection coefficient for an infinite resistive card, TE 
polarization. 
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Figure 4.38: Magnitude of the transmission coefficient for an infinite resistive card, 
TE polarization. 
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Frequency (GHz) 

Figure 4.39: Backscatter characteristics of the simple wedge geometry with impedance 
bulkhead, TM polarization. 

Figures 4.37 and 4.38 illustrate that as the angle of the incident TE wave ap- 
proaches grazing (0°), the transmission coefficient is large allowing energy to pass 
through the card. Thus, both tapered and constant resistive cards placed externally 
to the dielectric skin have little effect on the scattered field levels when the geometry 
is illuminated with TE polarization (see Figure 4.36). 

Impedance Bulkhead Surface 

Figures 4.39 and 4.40 illustrate the behavior of an impedance surface placed on the 
flat face of the planar bulkhead for TM and TE polarizations, respectively. The 
results indicate that continuing to increase the value of the impedance surface does 
not guarantee lower backscatter levels. 
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Figure 4.40: Backscatter characteristics of the simple wedge geometry with impedance 
bulkhead, TE polarization. 

4.4 Material Rod Nulling 

4.4.1 Dipole Arrays 

A demonstrated scattering reduction method is that of dipole array nulling. Dipole 
array nulling is effective in reducing scattering levels over a substantial bandwidth [12] 
[13]. Figure 4.41 ([12]) illustrates one example of the reduction possible by placing 
two co- linear dipole arrays in front of a leading edge model as described by Figures 
4.42 and 4.43. 

The goal of the following sections is to study the possibility of emulating the effects 
of lossy dipole arrays using perfect electric conducting or material rods. 

4.4.2 Perfect Electric Conducting Rods 

An illustration of the geometry(s) under consideration is shown in Figure 4.44. Rect- 
angular rods of various dimensions are placed in front of the conducting bulkhead as 
shown. There are an infinite number of combinations to explore when considering 
the number of rods, the rod dimensions, and the x and y location of the rods. As a 
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TOP VIEW 

Figure 4.42: Dipole array geometry. 
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REFERENCE PLOT GEOMETRY 

1 . rods 

a. radius of rods r roc j = .2 cm 

b. length of simulated ground plate dgp = 4 cm 

2. dipoles in front of rods 

a. inter-element separation D z = 5 cm 

b. length of elements 2L = 4.26 cm 

c. radius of elements r^jp = .025 cm 

d. resistance of elements R = 100 £2/cm 

e. offset dipole separation distance s = .2 cm 

f. dipole pair separation distance d = 1 cm 

g. ground plate spacing Dg = 2.25 cm 

Figure 4.43: Parameters of dipole array geometry. 

starting point, two rectangular rods with the same dimensions and distance from the 
conducting bulkhead are used. To limit the number of rod location possibilities (and 
computation time), the center of the rods are constrained to be located a maximum 
of 10 cm from the bulkhead in the x direction and within ±6.35 cm in the y direction 
(Figure 4.44). Additionally, the dimensions of the rods are set between 1 and 2 cm 
in length and height. 



Figure 4.44: Location and dimensions of material rods. 
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Frequency (GHz) 

Figure 4.45: Backscattered field from planar bulkhead with PEC rods placed in front 
of bulkhead, TM polarization. Solid - planar bulkhead alone, dashed - planar bulk- 
head with rods placed at x po ^=8 cm, y pom =±l cm, with A x =2 cm, and A„=l cm. 

The dimensions of the rods and each rods’ ±y location are held constant while 
varying the separation distance between the bulkhead and rods. The length, width, or 
±y location of the rods are changed for each scan of separation distance as a means of 
systematically investigating the effects of each parameter combination. The scattering 
signature of one of the more desirable responses is shown in Figure 4.45. It is possible 
that other combinations of parameters (i.e. a greater number or rods, larger rod size, 
larger separation distance, etc.) exist that are more effective in lowering the scattered 
field levels at the necessary frequencies. 

4.4.3 Lossy Dielectric Rods 

TM Polarization 

An initial consideration for this technique is the value of the dielectric constant used 
for the lossy dielectric rods. It is found that better nulling effects are obtained when 
the loss tangent of the material is large. An inspection of Table 4.1 shows that the 
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Figure 4.46: Backscattered field from planar bulkhead with four dielectric rods 

(e r =23.3-jl6.3) placed in front of bulkhead, TM polarization. Planar bulkhead alone 
- solid; planar bulkhead with rods placed at 3Cpo& n — f and 8 cm, ]/pojn”il cm, with 
A x =2 cm, A v =l cm - dashed. 

material under consideration with the largest loss tangent has a relative dielectric 
constant of e r =23.3-jl6.3. 

Using similar methods as described in the previous section to search for optimum 
rod placement and dimensions combined with results found in [13], a dielectric rod 
geometry is found that has desirable nulling effects. The scattering signature and 
dimensions of the geometry are given in Figure 4.46. Note that four dielectric rods 
are used to achieve the results. The first set of calculations are executed using two 
rods. Once the dual rod geometry with the best response is found, an additional two 
rods are added to investigate their effect. It is found that the two additional rods 
provide greater bandwidth. 

For the previous calculations, dielectric rods were considered with relative per- 
meability of fi r = 1. Figure 4.47 shows the effects of using the same configuration of 
dielectric rods as in Figure 4.46, but with /x r =4.8-j.4. This value of /i r is perhaps 
unrealistic, but it demonstrates that further reduction may be possible by altering 
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Frequency (GHz) 

Figure 4.47: Backscattered field from planar bulkhead with four dielectric rods 

(e r =23.3-jl6.3) placed in front of bulkhead, TM polarization. Planar bulkhead alone 
- solid; planar bulkhead with dielectric rods, /w«=l - short dashed; planar bulkhead 
with dielectric rods, fJ. ro d»= 4.8-j .4 - long dashed. 

the permeability of the rods. Note that a non-zero fi r is used with e r =23.3-jl6.3. It 
is found that varying fi r with e r =l is ineffective for this configuration. 

TE Polarization 

The desirable four rod configuration found in the previous section for TM polariza- 
tion did little to reduce backscatter levels when considering TE polarization. Many 
combinations of the four rod geometry were experimented with by varying the rod 
dimensions and placement, with negligible effect. Several two rod geometries were 
also considered, but still did little to influence the scattered field. 

Finally, a single material rod geometry was found to be effective. The resulting 
geometry is shown in Figure 4.48. The dimensions and location of the rod were 
found using values of e r =23.3-jl6.3 and fi T = 1. Figure 4.49 illustrates the backscatter 
signature of the geometry, as well as the improvement obtained by changing the 
relative permeability from n T =l to /z r =4.8-j.4. 
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Figure 4.48: Location and dimensions of material rod for TE polarization. 



Frequency (GHz) 


Figure 4.49: Backscattered field from planar bulkhead with single dielectric rod 

(e r =23.3-jl6.3) placed in front of bulkhead, TE polarization: Xpo tn =5 cm, A X =A V =8 
cm. Planar bulkhead alone - solid; planar bulkhead with dielectric rod, /i ro j= 1 - short 
dashed; planar bulkhead with dielectric rod, ^ ro j=4.8-j.4 - long dashed. 
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Chapter 5 

Composite Edge Termination 
Analysis 

5.1 Introduction 

Building upon the results of the previous chapter, the three planar bulkhead config- 
urations of Figure 5.1 are considered for further analysis. 

In the first configuration, internal material loadings are used in combination with 
resistive cards placed exterior to the dielectric skin. In the second configuration, 
the conducting bulkhead is shaped along with the use of resistive cards and internal 
loadings. The last configuration considers the nulling effects of material rods used 
in the presence of the dielectric skin and resistive cards. Both free space and lossy 
dielectric materials are considered for use in the region surrounding the material rods. 

5.2 Case I: Simple Wedge Geometry 

5.2.1 Resistive Cards Placed on Exterior Surface of Dielec- 
tric Skin 

The geometry studied in this section is described by Case I of Figure 5.1. The relative 
dielectric constant of the skin is set to e, kin =6, with free space in the internal region. 
The value of the resistive cards external to the dielectric skin are varied from 100 to 
500fi/D, with the signatures of the resulting geometries shown in Figures 5.2 and 5.3 
for TM and TE polarization, respectively. 


82 



Echo Width (dBm) 


45.2 cm 




Figure 5.1: Final geometries under consideration. 



Figure 5.2: Echo widths for simple wedge geometry with resistive card, TM polariza- 




Frequency (GHz) 

Figure 5.3: Echo widths for simple wedge geometry with resistive card, TE polariza- 
tion. 

Only the signatures for resistive card values of 100 and 500fl/D are included in the 
curves for both polarizations, since the variation between Q / □ values in this range 
is minimal. Resistive cards are more effective in reducing scattering levels for this 
configuration when illuminated with TM polarization. Additionally, the results of 
Section 4.3.3 help to illustrate this point. 

5.2.2 Constant Resistive Card Combined with Various Lossy 
Internal Loadings 

The previous section demonstrated that when considering TM polarization, the per- 
formance of a resistive card of 100Q/n had similar behavior as a resistive card of 
500n/D. For TE polarization, the difference in response between the two resistive 
card values is more dramatic, with the lOOfl/d card having slightly better perfor- 
mance. 

Based upon the results of the previous section, the geometry is altered by using 
a constant resistive card value of lOOfi/D while changing the internal region loading 
in accordance with the various e T values of Table 4.1. 
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Figure 5.4: Echo widths for simple wedge geometry with 10017/ □ resistive card with 
internal loading (e r ), TM polarization. 

Figures 5.4 and 5.5 give selected results for TM and TE polarization, respectively. 
For TM polarization, a resistive card placed on the outer surface of the dielectric skin 
limits the effectiveness of loading the interior region. This point is illustrated in Figure 
5.4 by providing the signature corresponding to an interior loading of e r =11.8-j3.16 
without the resistive card in place. For TE polarization, levels can be further reduced 
by loading the internal region. In addition, using an internal loading of e r =4.8-j.4 is 
comparatively as effective as an internal loading of e,.=11.8-j3.16. Thus, an internal 
loading of e r =4.8-j.4 will be used for calculations in the next section. Using a value 
of e r with smaller magnitude also reduces computation time as less unknowns are 
required to adequately describe the geometry. 

5.2.3 Constant Internal Loading Combined with Various 
Resistive Card Values 

Using the geometry of the previous section with an internal loading of e r =4.8-j.4 
and resistive card values in the range 100 - 500Q/D yields the backscatter signatures 
shown in Figures 5.6 and 5.7 for TM and TE polarization, respectively. There is 
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Figure 5.5: Echo widths for simple wedge geometry with lOOft/D resistive card with 
internal loading (e r ), TE polarization. 

little deviation in scattering levels between resistive card values in this range with 
TM polarization. This result again confirms that placement of a resistive card on 
the external surface of the dielectric skin reduces the effectiveness of the internal 
loading. For TE polarization, a resistive card value of 200Q/D is most effective in the 
frequency range from .5 to 1 GHz, with a value of lOOfl/D becoming most effective 
in the range 1 to 3 GHz. 

5.3 Case II: Simple Wedge Geometry with Shap- 
ing of Conducting Bulkhead 

5.3.1 Resistive Cards Placed on Exterior Surface of Dielec- 
tric Skin 

The geometry examined in this section is described by Case II of Figure 5.1. The 
difference from Case I is in the shape of the conducting bulkhead. The front portion 
of the bulkhead is extended to form a wedge as shown, with the distance from the 
bulkhead junction set at 25.4 cm. Once again, the relative dielectric constant of 
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Figure 5.6: Echo widths for simple wedge geometry with resistive card, internal load 
ing e r =4.8-j.4, TM polarization. 



Figure 5.7: Echo widths for simple wedge geometry with resistive card, internal load 
ing e r =4.8-j.4, TE polarization. 




Figure 5.8: Echo widths for shaped conducting bulkhead geometry with resistive card, 
TM polarization. 

the skin is set to e < fc tn =6 ) with free space in the internal region. The value of the 
resistive cards external to the dielectric skin are varied from 100 to 500fi/d, with the 
signatures of the resulting geometries shown in Figures 5.8 and 5.9 for TM and TE 
polarization, respectively. 

For TE polarization, resistive cards in the specified range are largely ineffective 
when used in conjunction with a wedge shaped conducting bulkhead. For TM po- 
larization, resistive cards help reduce scattering levels at low frequencies, with only 
slight signature variation between card values. 

5.3.2 Constant Resistive Card Combined with Various Lossy 
Internal Loadings 

The results of the previous section demonstrate that there is little difference in scat- 
tering signatures between resistive card values. The calculations performed in this 
section utilize a resistive card of lOOfi/D while varying the internal region loading. 
The results for TM and TE polarizations are given in Figures 5.10 and 5.11, respec- 
tively. 
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Figure 5.9: Echo widths for shaped conducting bulkhead geometry with resistive card, 
TE polarization. 



Figure 5.10: Echo widths for shaped conducting bulkhead geometry with 100fi/O 
resistive card, TM polarization. 
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Figure 5.11: Echo widths for shaped conducting bulkhead geometry with 100Q/D 
resistive card, TE polarization. 

As was found for Case I, a resistive card placed on the outer surface of the dielec- 
tric skin limits the effectiveness of loading the interior region when considering TM 
polarization. There is, however, a small deviation between the TE signatures of the 
bulkheads in Cases I and II when considering a constant resistive card of lOOft/D 
with various internal region loadings. 

In Case I (TE), loading the interior region is effective in further influencing the 
backscattered fields when used in conjunction with a resistive card on the external 
surface of the dielectric skin. In contrast, loading the interior region of the geometry 
of Case II (TE) does not have as dramatic of an effect. A possible explanation for 
this behavior is that the wedge shaped bulkhead alone already provides significant 
reduction when compared to the signature of the planar bulkhead alone. 

5.3.3 Impenetrable Impedance Bulkhead 

Section 4.3.3 investigates the effects of placing an impedance surface on the flat por- 
tion of the planar conducting bulkhead. Applying the same concept to the wedge 
shaped portion of the conducting bulkhead and considering only the bulkhead in free 
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Figure 5.12: Echo widths for wedge shaped bulkhead with impenetrable impedance 
surface on wedge portion of bulkhead, TM polarization. 

space yields the results shown in Figures 5.12 and 5.13 for TM and TE polarization, 
respectively. The bulkhead is considered in absence of the dielectric skin to illustrate 
that it is the shape of the bulkhead that limits the effectiveness of the impedance 
surface for both polarizations. 

5.4 Case III: Simple Wedge Geometry with Place- 
ment of Dielectric Rods 

5.4.1 Dielectric Rod Placement Analysis 

This section builds upon the dielectric rod geometry of Figure 4.46 in Section 4.4.3. 
The difference between the geometry of Section 4.4.3 and the current geometry is 
that the bulkhead and rods are now considered in the presence of the dielectric skin, 
as shown in Figure 5.15. 

In deciding what material should be used for the rods, swept frequency calcu- 
lations are run for rod relative dielectric constants of Table 4.1. The results are 
shown in Figure 5.14 for TM polarization. Also included is a calculation using PEC 
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Frequency (GHz) 


Figure 5.13: Echo widths for wedge shaped bulkhead with impenetrable impedance 
surface on wedge portion of bulkhead, TE polarization. 

rods in place of lossy dielectric rods. The results illustrate that using dielectric rods 
with a high loss tangent (e r =23.3-jl6.3) is more effective in reducing low frequency 
backscatter than using PEC rods. 

The GA is applied to search for an optimum combination of XI, X2, and Y, ep (as 
shown in Figure 5.15) with discretizations given by the GA template file of Figure 
5.16. The search is performed using three different fitness functionals corresponding 
to TM echo width at .5, .7, or 1 GHz. 

The best results for the three GA runs are summarized in Tables 5.1 through 5.3. 
The signature for each of the three search cases have varying performance as shown 
in Figures 5.17 through 5.19. 

As a means of lowering the TE return without largely affecting the TM return, a 
500fi impedance is applied to the flat portion of the conducting bulkhead with the 
resulting signatures shown in Figures 5.20 through 5.22. Unfortunately, as was noted 
in Section 4.3.2, this impedance surface may not be realistic. 
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genes : 3 

gene 0 
min : 2 
max : 5.75 
values: 16 
format : '/, 9 . 3f 

gene 1 
min: 3 
max: 6.75 
values: 16 
format: '/,9.3f 

gene 2 
min: 1 
max: 4.75 
values: 16 
format : ’/,9 . 3f 


Figure 5.16: “template” file corresponding to the search for optimum dielectric rod 
placement. 
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Table 5.1: Search results for optimum combination of XI, X2, and Y sep with Echo 
Width evaluated at .5 GHz, TM polarization. 




Y 

1 aep 

Echo Width 

Generation 

trial 



(cm) 

(dBm) 





1.000 

-1.3004e+01 

6 

168 

5.750 

3.000 

1.250 

-1.2588e-f 01 

15 

262 

5.250 



-1.2926e+01 

2 

79 

4.750 



-1.2827e+01 

8 

193 

IHH 



-1.2912e+01 

9 

205 

5.750 



-1.2954e+01 

10 

221 

5.500 

3.000 

1.000 


3 

106 



1.000 

-1.3510e+01 

2 

84 



1.000 

-1.2637e+01 

11 

239 

5.000 

3.000 

1.000 

-1.2730e+01 

1 

61 
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Table 5.2: Search results for optimum combination of XI, X2, and Y tep with Echo 
Width evaluated at .7 GHz, TM polarization. 


XI 

(cm) 

X2 

(cm) 

Y 

J eep 

(cm) 

Echo Width 
(dBm) 

Generation 

trial 

2.500 

3.250 

1.250 

-3.2803e+01 

0 

17 

2.500 

3.000 

1.500 

-3.0087e+01 

5 

147 

2.250 

3.500 

1.500 

-3.4257e+01 

13 

220 

2.750 

3.000 

1.500 

-3.8261e+01 

12 

211 



1.250 

-3.8599e+01 

3 


2.750 

3.250 

1.500 

-3.0989e+01 

0 

43 

2.750 

3.250 

1.250 

-3.2524e+01 

4 

144 

2.250 

3.500 

1.250 

-3.1964e+01 

12 

209 

2.500 

3.500 

1.000 

-3.0920e+01 

2 

94 

2.500 

3.250 

1.500 

-3.4025e+01 

1 

62 
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Table 5.3: Search results for optimum combination of XI, X2, and Y, ep with Echo 
Width evaluated at 1 GHz, TM polarization. 


XI 

X2 

Y 

1 aep 

Echo Width 

Generation 

trial 

(cm) 

(cm) 

(cm) 

(dBm) 



4.750 

3.500 

1.000 

-1.8635e+01 

10 

241 

4.500 

3.750 

1.000 

-1.9272e+01 

9 

238 

4.750 

3.250 

1.000 

-1.8313e+01 

3 

109 

3.500 

4.000 

1.000 

-1.8340e+01 

20 

400 

4.250 

3.750 

1.000 

- 1 * 8898e —(- 0 1 

12 

281 

3.750 

4.000 

1.000 

-1.8575e+01 

17 

358 

4.500 

3.500 

1.000 

-1.9361e+01 

12 

286 

4.000 

3.750 

1.000 

-1.8373e+01 

5 

142 

4.250 

4.000 

1.000 

-1.8442e+01 

20 

387 

4.000 

4.000 

1.000 

-1.8476e+01 

4 

114 


97 





Echo Width (dBm) 3 Echo Width (dBm) 



Frequency (GHz) 

5.17: Echo Widths for dielectric rod geometry optimized to 0.5 GHz. 



Frequency (GHz) 

Figure 5.18: Echo Widths for dielectric rod geometry optimized to 0.7 GHz. 


98 




Echo Width (dBm) Echo Width (dBm) 



Frequency (GHz) 


5.19: Echo Widths for dielectric rod geometry optimized to 1.0 GHz. 



Frequency (GHz) 


Figure 5.20: Echo Widths for dielectric rod geometry optimized to 0.5 GHz, with a 
500H impedance placed on the flat portion of the conducting bulkhead. 





-60 3 

0 0.5 1 1.5 2 2.5 3 

Frequency (GHz) 


Figure 5.21: Echo Widths for dielectric rod geometry optimized to 0.7 GHz, with a 
500fl impedance placed on the flat portion of the conducting bulkhead. 
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Figure 5.22: Echo Widths for dielectric rod geometry optimized to 1.0 GHz, with a 
500fi impedance placed on the flat portion of the conducting bulkhead. 
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5.4.2 Variation of Surrounding Dielectric Material 

This section illustrates the effects of loading the (internal) region surrounding the 
dielectric rods with various lossy dielectric materials of Table 4.1. Figures 5.23 and 
5.24 give the results for TM and TE polarization, respectively. The backscatter 
characteristics of the resulting geometry are similar to the performance obtained when 
loading the internal region with a homogeneous lossy dielectric material as described 
in Section 4.3.2. 
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Chapter 6 
Conclusions 


In this report, a variety of edge termination methods were studied. It was found that 
the planar bulkhead with skin geometry has the most desirable backscatter character- 
istics for both polarizations provided the internal region is loaded with a sufficiently 
lossy dielectric material as shown in Figures 4.31 and 4.32. Placement of material 
rods also has desirable backscatter reduction characteristics for both polarizations 
if a surface impedance can be placed on the flat portion of the planar conducting 
bulkhead, as shown in Figures 5.20 through 5.22. 

A method that did not work for both polarizations is the use of resistive cards. 
Resistive cards on the external surface of the dielectric skin were found to be an 
undesirable scattering reduction method since the presence of the card limits the 
effectiveness of loading the internal region for TM incidence as shown in Figure 5.4. 
Also, shaping of the internal conducting bulkhead reduced (TE) scattering from the 
tip and bulkhead junction, but did little to lower the overall backscatter levels of the 
bulkhead. 
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